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Abstract 



The thermalization time for a Quark-Gluon-Plasma is estimated from a quan- 
tum transport model beyond the quasi-particle approach (or kinetic gas the- 
ory). While our ansatz is crude concerning the properties of "real" quarks 
and gluons, it nevertheless takes very serious the basic principles of quantum 
field theory for non-equilibrium states. It is found, that the thermalization 
time obtained from quantum transport theory is substantially longer than 
from kinetic theory. In our view this casts some doubts on scenarios which a 
priori assume a thermalized quark gluon plasma. 
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I. INTRODUCTION 



In the past decade the effort invested in ultrarelativistic heavy ion collisions (URHIC) has 
grown considerably |TJ . The general hope is, that at some time in the near future one may be 
able to observe an excursion of strongly interacting matter from the state of hadrons before 
the collision into the phase of a quark-gluon plasma (QGP). Consequently, the discussion 
of possible signals from such a shortlived state is quite vivid: Weakly interacting probes like 
photons or lepton pairs, as well as strongly interacting signals like those presented by quark 
flavors of higher mass have been proposed. Similar to most of these investigations is the 
assumption of a thermalized plasma phase, followed by the calculation of the time evolution 
along one or the other line of physical reasoning. 

With the present paper we address the question, whether such a thermalized phase is 
reached at all. We limit ourselves to the physical scenario that one may reach in future 
URHIC: A sea of gluons, initially at low temperature, is heated to a very high temperature 
over a short time. In this hot glue, quark-antiquark pairs are popping up - until at the very 
end a thermal equilibrium in the sense of a degenerate plasma is reached. Of great interest 
also for the future experiments is the time scale for this equilibration. 

If one considers quarks and gluons in a standard (Boltzmann-like) transport theory, this 
time scale is obtained as ~ 1 fm/c - short enough to reach a pseudo-equilibrium in the time 
span available in heavy ion collisions. However, we have serious doubts that the requirements 
for the applicability of kinetic gas theory are fulfilled in a QGP: The thermal scattering of 
constituents occurs so frequently, that subsequent collisions overlap quantum mechanically. 

In other words, we expect that in a QGP the off-shell propagation of particles plays 
an important role and therefore the quasi-particle approximation is not applicable. For- 
mally, we describe this by including a nontrivial spectral function of the constituents, i.e., 
in the description we go beyond the quasi-particle approximation. However, the spectral 
broadening of particles in a hot plasma, although described on the same footing as a gas of 
resonances, does not imply that our particles may decay. Their continuous mass spectrum 
merely represents the thermal scattering by the other components of the system. 

For the limited purpose of the present paper, we furthermore make some physically 
motivated assumptions: 

1. We assume, that the self energy function for the quarks is dominated by gluonic 
contributions. This is justified because the quark-quark scattering cross section is 
much smaller than the quark-gluon cross section. 

2. The gluon background is dominated by external conditions, i.e., we neglect the back- 
reaction of quarks on the gluon distribution. 

3. The external conditions determining the gluon field are changing in a short time in- 
terval, and the system is translationally invariant in 3-dimensional coordinate space. 

For these assumptions, we also have a practical reason: They allow for a clean separation 
of various aspects of the quantum transport problem, whereas this separation is difficult 
(if not impossible) when considering more realistic systems. While our primary motivation 
to go beyond the quasi-particle picture therefore is a physical one, we may also adopt a 
mathematically rigorous stance. The Narnhofer-Thirring theorem || states, that interacting 
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systems at finite temperature cannot be described by particles with a sharp dispersion law. 
Ignoring this mathematical fact one finds as an echo serious infrared divergences in high 
temperature quantum chromodynamics (QCD). Consequently, these unphysical singularities 
are naturally removed within the approach of finite temperature field theory with continuous 
mass spectrum 0. 

The effect of off-shell propagation on relaxation processes was first investigated by 
Danielewicz M. He found a substantial slowdown of the equilibration, when solving the 
full quantum transport equation as compared to the solution of its kinetic approximation. 
In later years several attempts were made to unify transport theory with concepts going be- 
yond the quasi-particle picture The connection between transport theory and quantum 
field theory with continuous mass spectrum was established in refs. |7||| . A first application 
of these principles to the QGP relaxation problem P,|T0fl points into the same direction as 
the results of Danielewicz. 

The paper is organized as follows: In the next section we give a brief introduction into 
the formalism necessary for non-equilibrium quantum fields. In section III we discuss our 
approximate spectral function, followed by a solution of the quantum transport equation 
in section IV. In section V we show, how one may derive a transport equation beyond the 
Boltzmann-like quasi-particle approximation from the equations in section III, and we solve 
it in section VI. Conclusions are drawn in the final section of the present work. 



II. MATRIX- VALUED SCHWINGER-DYSON EQUATION 

As has been pointed out by various authors, the description of dynamical (time depen- 
dent) quantum phenomena in a statistical ensemble necessitates a formalism with a doubled 
Hilbert space [f|,[7||- For our purpose the relevant content of this formalism is that its two- 
point Green functions are 2x2 matrix- valued. We leave it to the reader to chose either 
the conventional Schwinger-Keldysh, or Closed-Time Path (CTP) Green function formalism 
11,12], or the technically simpler method of Thermo Field Dynamics (TFD) [|1^ . 



Within this matrix formulation, we consider the Schwinger-Dyson equation for the full 
quark propagator 

S = S + S S S . (1) 

Here So is the free and S the full two-point Green function of the quark field, E is the full 
self energy and the generalized product of these is to be understood as a matrix product 
(thermal and spinor indices) and an integration (each of the matrices is a function of two 
space coordinates): 

s S/0^ = E \ d 'v^% s t- ( 2 ) 

3 

Throughout this paper we use the convention to write space-time and momentum variables 
also as lower indices, e.g. Ti xy = E(x, y). 

For the purpose of treating transport equations, we switch to the mixed (or Wigner) 
representation of functions depending on two space-time coordinates: 
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t XP = jd\x - y) exp {iP^x - yY) E xy . (3) 

The ~-sign will be dropped henceforth. The Wigner transform of the convolution £ G as 
defined in (||) is a nontrivial step: Formally it may be expressed as a gradient expansion 

jd\x - y) exp - y)") = t XP exp (-i G X p • (4) 

O is a 2nd order differential operator acting to the left and to the right side, 

3= ^ (d x dp - dpd x ) . (5) 

The arrows on the derivatives display the directions in which they act. Explicitly, the 
first-order term is 

v ~ i (dY, X pdG X p dY, X pdG XP \ 

Sxp OG XP = -\^ - d ^ j , (6) 

i.e., this is the Poisson bracket of the two quantities G and E. 

In the CTP formulation as well as in the a = 1 parameterization of TFD [H|, the matrix 
elements of S, the time-ordered propagator S 11 , the anti-time-ordered propagator S 22 and 
the IP-correlations S 12 and S 21 , are not linearly independent 

S 11 + S 22 - S 12 -S 21 =0. (7) 

A similar relation holds for the free propagator So and with different signs for the components 
of the self energy: 

S ll + ^22 + S 12 + E 21 = (g) 

These relations lead to different representations of the retarded and advanced propagators 
and self-energies 

qR _ nil q12 _ q21 c22 qA _ oil o21 _ ol2 r-22 

E R = S 11 + S 12 = -S 21 - S 22 , S A = S 11 + S 21 = -S 12 - S 22 . (9) 

Due to the above linear relations, the four components of the Schwinger-Dyson equation are 
not independent, the matrix equation can be simplified by a linear transformation. This lin- 



ear transformation, which one may conveniently express as a matrix transformation fl5| , |14 
has a physical interpretation only in the TFD formalism, see ref. 0. The transformation 
matrices B are 

7). w»>>MLi (!-„))• (10) 

depending on one parameter only. Due to the linear relations (0) and (|8|), we obtain for any 
value of this parameter 
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with r 3 = diag(l, —1). 

For equilibrium systems the off-diagonal elements of the original matrix valued propaga- 



tor fulfill the Kubo-Martin-Schwinger (KMS) condition JIB 

(1 - n F ( Po )) S 12 (p ,p) + n F (p )S 21 (p ,p) = , (12) 



where n F (E) is the Fermi-Dirac equilibrium distribution function at temperature T, 

1 



ME) = ^ , ; • (13) 



It is obvious, that therefore in equilibrium systems the transformation (|Tl|) with n = n F 
results in complete diagonalization of the matrix S [I4|,7|. The same distribution function 
may also be used to diagonalize the matrix valued self energy function, because it fulfills a 
relation similar to flT2l). 

In the next step the above condition is generalized to 

(1 - N XP ) S 12 P + N XP Sf P = , (14) 

which defines the non-equilibrium distribution function N X p depending on coordinates and 
momenta. In other words, we choose as non-equilibrium distribution function the parameter 
which diagonalizes the non-equilibrium propagator under the transformation (|TT|) : 



B{N XP ) r 3 S XP {B{N XP )Y l = ( S * p ^ ) . (15) 

Generally however, this parameter N XP does not diagonalize the non-equilibrium self en- 
ergy function. For this purpose one has to introduce another parameter N^ p , and the 
transformation of the matrix valued self energy function then yields 

B(N XP ) S XP r 3 (B(N XP )r = ( ^ P 2 ^ N xp-Nxp)T xp ^ (lg) 

where T is the (common) imaginary part of retarded and advanced self energy. From this 
splitting follows that the off-diagonal element in Eq. vanishes for equilibrium states, 

because then N XP = N XP = np(po). This we may use as an indication that the off-diagonal 
term of ( |TED is the "collision term" of a non-equilibrium system. 

The splitting into real and imaginary part is useful for self energy as well as propagator, 

S XP = G XP =f inA XP 

E x 'p = ReS XP =f inT XP . (17) 

Each of these functions is assumed to be real, and they are spinor-valued. The part of the 
self energy function which is local in space and time, is in our notation included in ReS^p- 
For a relativistic model this is the Hartree part of S. 
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From eqns. ( |i5|) and (|16|) follows, that the diagonal elements of the matrix transformed 
Schwinger-Dyson equation are retarded and advanced Schwinger-Dyson equation. The off- 
diagonal element is a transport equation, because it determines the non-equilibrium distri- 
bution function. As said above, the driving force of the transport part is the deviation of 
the non-equilibrium distribution function from the corresponding equilibrium distribution 
function to which the whole system tends. 

A special virtue of this formulation is the fact, that with the KMS condition (|T2p we 
have already shown that this non-equilibrium distribution function reaches a Fermi-Dirac 
distribution in case we are considering an equilibrium system. 

We now consider the equations obtained by action of Dirac differential operators (= 
inverse free propagators) on the matrix-transformed Schwinger-Dyson equation. They act 
on the free retarded and advanced propagator as 

So 1 S$ P = (V 7 , -M+ ~ 7/ ^) S$ P = 1 . (18) 
Acting with the inverse free propagator on the full Schwinger-Dyson equation, we obtain 

So 1 S A X ? = 1 + E$? exp (-i S) . (19) 

Adding and subtracting the equations for the retarded and advanced propagators, we obtain 
equations of motion for the spectral function and the real part of the retarded and advanced 
propagator 

2ivr5 - 1 = T, A XP exp (-i 3) S XP - exp (-i o) S$ P , (20) 

G x 'p = 1 + ± (E£ p exp (-i S) Sl P + Z XP exp (-i o) S XP ) (21) 

for the diagonal components of the transformed Schwinger-Dyson equation. The solution in 
zero order O leads to the solution for a equilibrium system. The equations in first order of 
the gradient approxiamtion give a transport equation for the spectral function which is in 
equilibrium solved by the solution in zero order O. 

A very important fact about these equations has to be emphasized ||: In general, they 
do not admit a ^-function solution for the spectral function A X p even in zero order of the 
gradient expansion. This has led to much confusion in papers deriving transport equations 
from the Schwinger-Dyson equation. We therefore state very clearly what is obvious for 
mathematical as well as physical reasons: There is not such thing as a mass shell constraint 
in quantum transport theory ! To insert such a constraint into any model by hand must be 
considered dangerous unless it is tested against the full calculation. 

The off-diagonal component of the transformed Schwinger-Dyson equation reads, after 
acting on it with the inverse free propagator, 

S^Sftp = X£ P exp (-i O) Sfp - Sf P exp (-i S) S XP (22) 

(see 0, pp.307] for details). A similar equation holds for the propagator S 12 

S^S 1 ^ = p exp (-i 5) S 1 ^ + T}' P exp (-i o) S XP (23) 
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Using the equation for the spectral function (eq. fl20|)) and S%p = Zi^NxpAxp, we obtain 
a differential equation for Nxp 

- Tr^Axpd x N XP = P exp (-i $) + S^ 2 P exp (-i $) 

+iV X p(s5 P exp (-i O) S-fp - S^p exp (-i o) s£ p ) . (24) 

The transport equations are usually derived by the terms up to first order O, but it is visible 
that transport equations contain in principle arbitrary orders of the operator O. Below for a 
certain approximation the transport equation will be solved up to all orders in the gradient 
expansion. 

III. ANSATZ FOR SELF ENERGY AND SPECTRAL FUNCTIONS 

As we have seen above, the matrix valued propagator has only three independent com- 
ponents, two of which are furthermore complex conjugate. We have also shown how one 
may use the KMS condition to eliminate the off-diagonal component of this propagator in 
favor of the distribution function parameter. 

We now tackle the remaining two pieces, i.e., we deal with retarded and advanced prop- 
agator. In coordinate space they are zero in the forward resp. backward lightcone (outside, 
both are identically zero), 

S R y A = T 27ri9 (±0c o - y )) A xy . (25) 

Consequently, their Wigner representations are analytical functions in the upper resp. lower 
complex energy halfplane. Hence, even for non-equilibrium states we may write in the mixed 
(or Wigner) representation 

oo 

S R > A (E, p, X) = G XP T TTiAxr = J dE' f^'^fl ■ (26) 

— oo 

Without elaboration at this point we emphasize, that one may not calculate the real part 
of the quark propagator by a method different from eq. (p6|), e.g. perturbatively. In a 
relativistic system, doing so may lead to unphysical poles in the retarded propagator, on the 
wrong side of the real energy axis - the retarded propagator then does not fulfill eq. (p5|). 

The principal difference between equilibrium and non-equilibrium states at this point are 
the properties of the function 

A(E,p) = T -MS R > A (E,p)) = -L (S A {E,p) - S R (E,p)) . (27) 

7T Z7T1 

In equilibrium states it has spectral properties, i.e., it is normalized and positive semidefinite. 
In non-equilibrium states positivity is not guaranteed, we can only use the term "spectral 
function" in a generalized fashion. Normalization of the 7°-component of A according to 

oo 

JdETi[-f°A(E,p)]=2 (28) 
o 
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is guaranteed also in non-equilibrium states: It is a consequence of the canonical anti- 
commutation relations for the fields. 

For equilibrium states we may now combine the facts into a very compact notation for 
the matrix valued propagator, using the Bogoliubov matrix defined in ([K]) fl4]j7| 



S (ab) (Po, P )= JdEr,{B{n F {E)))- 1 ho~E + ie 



B{n F {E)) . (29) 



/ 

Po — E — le / 



For free fermions of mass m the spectral function is 

A(E, p) — ► (£7° + P7 + m) sign(E) 5(E 2 - p 2 - m 2 ) . (30) 

Instead of attempting a fully self-consistent numerical calculation of the spectral function 
for non-equilibrium states, we make another approximation. 

It is assumed, that the spectral function of the interacting system does not differ too 
much from the spectral function of quasi-particles. In particular, we make the ansatz: 

7t 7° (E 2 + uJ 2 t +l 2 t ) + 2Ejp + 2Em t 



(£2 - a;? - 7 f ) 2 + 4£? 7 f 
1 / co> t 7° + p7 + m t -tutl + PI + m t 



A{E,p,t) 



A-Kiut \ E — u t — yy t E + uo t - vy t 

ujg + p~f + m t +p~f + m t \ 

E - uj t + vy t E + u t + ht J 

where mt is the effective mass mt = m — ReS(p = 0), see below. Hence, we approximate the 
quark spectral function by two time-dependent parameters u t and 7t, which we may interpret 
as effective quasi-particle energy and effective spectral width and describe the energy pole 
as a function of time. One may argue about the validity of this approach, in particular 
whether not a momentum dependent spectral width is an absolute necessity for a realistic 
calculation. 

Let us discuss this in three steps. First of all, in ref. |17| the same approximation was 
used to obtain numerical results for photon radiation rates from a hot QGP. It was found, 
that in the proper temperature and energy regions these results agree very well with the 
hard thermal loop approximation scheme of QCD. 

Secondly, we may safely assume that the quarks appearing in the hot medium are slow 
- hence the properties of the quark distribution may be safely approximated by those of 
quarks at rest. Third, and more important from the aspect of model consistency is the fact 
that the expectation value of the anti-commutator of two quark fields is 

oo 

{^(x)My)}) = JdE l0- 3 A[E,p,^±^)exp(~i(E(x°-y o )-p(x-y))) . 

— oo 

(32) 

With our ansatz for the spectral function it is easy to show, that this commutator vanishes 
for spacelike distances x — y also in case Ut and jt are time-dependent [18| . Conversely, for 
a general momentum dependence of the parameters u and 7 this cannot be guaranteed. 
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By inspection of eqs. fl20|) and (|4|) we find, that beyond this only a quark self energy 
function is needed for a full determination of the function Axp- This self energy function 
is in general a functional of Axp again - which then leads to a complicated set of integro- 
differential equations for the self consistent determination of the retarded and advanced 
propagator. However, for the limited purpose of the present paper we have formulated 
simple approximations in the introductory part of the present paper. 

In our equations these ideas are introduced through replacing the X, P dependence of 
the self energy completely by a dependence on X = t, 

ReS Xo = ReS f = ReE f 9(t) + ReS^(-t) , 

T Xo =T t = g{T)>f(TfO{t)+T0(-t)) . (33) 

In other terms, we take as the physically most important pieces of the self energy function 
the Dirac scalar real part, and the Dirac "vector" imaginary part. The parameters ReE,, 
Re£/, denote the constant quark mass shift before and after time t — 0, Tj and Tf denote 
the corresponding temperatures. g{T) is some numerical factor given below. 

This ansatz corresponds to instantaneous heating of the gluon background and describes 
quarks in the corresponding complex external field. It is different from a previous treatment 
off the relaxation problem [IIJ, where only the imaginary part of the self energy function 
was used as input. 



IV. SOLUTION OF THE RETARDED AND ADVANCED EQUATION 



In the following we discuss the solutions for the time-dependent functions uj t and 7t from 
our ansatz for the spectral function, which are obtained when inserting the ansatz for the 
self-energy. Indeed, using this self energy it is possible to solve the r.h.s. of the equation of 
motions of the retarded and advanced functions to all orders in the gradient expansion, see 
appendix. In particular, the equations obtained after taking the trace over the Dirac indices 
and adding the real parts of Eq. (|20|) resp. subtracting the imaginary parts are 



1 



-ReTr 



-ImTr 

4 



j^-M+'-j^) (S R P + S XP ) 



+ -ReTr 
4 



xp exp -l 



°XP "i ^ 



xp exp -i 



7/i P" - M + i 7M 3H {Sxp-Sxp) 



-ImTr 

4 



xp exp 



-i O ] S R P - S£ P exp l -i 



qA 



(34) 



The complete expressions are given in the appendix. The choice of the above combinations 
of taking the real or imaginary part and adding or subtracting the equations for S R and 
S A ensures that drift term contributions for S R and S A resulting from the derivative on 
the l.h.s. drop out. Thus, only the zeroth order terms in the gradient expansion remain on 
the l.h.s.. Taking the unphysical limit of vanishing spectral width, the above equations lead 
to the mass shell constraint. The other two, possible combinations pick out the drift term 
contributions and lead to transport equations for the real and imaginary part of S R and S A . 
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For a system on non-relativistic particles these transport equations are fulfilled automatically 
in gradient expansion if one inserts the solution for S R and S A from the equations of zeroth 
order in the gradient expansion |7|,[T3|. Here, taking into account all orders of the gradient 
expansion, these equations lead to two additional equations which must be fulfilled and are 
discussed below. 

Performing further reductions, one obtains a set of two coupled nonlinear equations for 
7t and Up For a concrete evaluation, the energy parameter is chosen as E = uop This implies 
that the energy parameter is equal to the real part of the quasi-particle pole, the equations 
are thus evaluated at the peak of the spectral function. The result is: 



CO, 



2 2 u tlt (E 1 + F 1 )-(2u 2 + 1 2 t )(E 2 + F 2 ) 
p +m t H 



It = TrgT 



2o^(Si + Fi) - lt (E 2 + F 2 ) 



(35) 



with functions Ei >2 and F 12 defined as 
m 4 (ReSj — ReSj 



E, = - 



4u* 



6{t)e- 2ltt {2u;a t cos(4^t) + 7 2 sin(4j t i)) 



+6(-t)e 2 " ltt (2uja t cos(Auj t t) - 7 2 sin(Mt)) 



4u t 



e{t)e- 2ltt {Auj 2 + 7 2 - 7 2 aaa{Au t t) + 2u^ t sin(4^t)) 



+6(-t)e 2ltt (4uj 2 + 7 2 - 7 t 2 cos(4u; t t) - 2u tlt wn(4w t t)) 



9{t)e~ 2lt \Auj 2 + 7 2 + 7 2 cos(4o; t t) - 2u^ t sin(4o; t t)) 



+6(-t)e 2 " /tt (4uj 2 + 7 2 + 7 2 cos(Auj t t) + 2u t ~/ t sm(4u t t)) 



4 



- 9{t)e~ 2ltt {2uja t cos(4^t) + 7 2 sin(4c^)) 
+6(-t)e 2ltt {2uja t cos(4o; t t) - 7 2 sin(4^t)) . 

First we discuss limiting cases of the above equations. For t = ±oo we obtain 



(36) 



0J f 



=±0O = p 2 + (m + ReE /A ) 2 = p 2 + mJ /4 



7t=±c 



(37) 



Consequently, our equations have the proper boundary conditions for the effective quark 
mass and the spectral width of the quarks. 

Furthermore, one sees that at t — the effective width remains smooth and takes the 
value 7t=o = ng(Tf — Tj)/2, whereas the function uu t has a jump at t = 0: The limit from 
the left is uj 2 =0 _ = p 2 + mj(mj + m/)/2 and from the right is w 2 =0+ = p 2 + m/(mj + m/)/2. 

In each of these equations it is easy to distinguish between the contribution of quarks 
and antiquarks: If one takes into account only the positive energy pieces of the spectral 
function, terms oscillating with frequency 4u t drop out. 
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We now turn to a numerical evaluation of these formal results. To this end, we first have 



to specify what we mean by the factor g in eq. (33). In our view its is not clear to the 
present date, how this parameter is related to the strong coupling constant: Some authors 
use a quadratic, others a linear dependence on the actual g$. A thorough discussion is 



carried out in the second paper of ref. JI7|]. Since it was obtained there, that the numerical 
results are quite independent of the actual functional structure, we adopt a pragmatic view 
in the following and use a simple expression estimated within the hard thermal loop scheme 
assuming a two flavor system and quarks with zero momentum (p 2 = 0) p0| : 

5 63 

7(T/)= 9 2 T f = 1.88 a s T f , 
g = 0.598 a s ■ (38) 

The temperature dependent coupling g follows a parametrization of the strong coupling 
constant as = (? 2 /47r given by Karsch ||21||: 



as(T) = (33 - 2n / 6 ) 7 ln(8T/T c ) ' (39) 

where T c is the critical temperature and rif the number of flavors. Some values for the 
coupling g with the corresponding temperature and critical temperature are listed in Table 
1. For temperatures less than T c the spectral width of the quarks is chosen infinitely small, 
in agreement with more elaborate self-consistent calculation within a generalized Nambu - 
Jona-Lasinio model [|17j . 



The effective quark mass for temperatures less than T c is taken as the standard con- 
stituent quark mass mj=300 MeV. Above the critical temperature, we take a termal quark 
mass from |E| 

2 T 2 2 2 

m T = —T~ + ITT ■ 40 

The term containing the chemical potential is small in comparison with the first term and 
neglected in the calculations. The thermal masses are shown in Table 2. 

Plugging these parametrizations into our highly nonlinear equations then requires some 
effort to solve them numerically. Let us now comment on the results of this model calculation. 
In Fig. 1 m t = u t (p = 0) is depicted as a function of time. Furthermore for the special case 
of Tj = T c the temperature dependence of the solution is indicated. As shown above, the 
constitutent quark mass of 300 MeV is obtained as the boundary value. Approaching t = 
the effective quark mass decreases with time until it reaches the left-sided limiting value at 
t = 0_ discussed above. Then is jumps down to the right-sided limiting value at t = + and 
decreases further and approaches the final value for the quark-gluon plasma. 

Conversely, the spectral width gamma increases with time, from the small value in the 
nuclear environment to a large value in the QGP. Its value is shown in Fig 2. Note, that due 
to our ansatz functions we find a jump in the effective quark mass, whereas jt is smooth at 
t = 0. Both curves exhibit quantum oscillations with the frequency AuJt, becoming slower 
due to the decrease of u t as a function of time. 

It is furthermore necessary to point out that the effects for t < are no acausal behaviour. 
They arise in the present formulation from the fact that t is the center-of-mass time of a 
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2-point function and not the physical time of the system. Therefore, even for t < there 
are contributions coming from the time after the jump in temperature since it is possible to 
choose an adapted relative time. Naturally these contributions decrease to zero if t — > — oo 
as it can be seen in Fig 1 and 2. 

Besides the above two equations we can get two further equations for the spectral function 
by adding the imaginary parts and by subtracting the real parts of the equations for S R and 
S A . These equations could in principle be solved if we would not have fixed the shape of the 
spectral function as function of E and p by the above ansatz. Consequently, in our picture 
they are fulfilled only approximately. 



V. TRANSPORT EQUATION 



We now turn to the primary task of the present paper, i.e., to the evolution equation for 
the quark distribution function. This transport equation is obtained from eq. @ by taking 
the trace over spinor indices and separation into real and imaginary part: 



ReTr\ - 7r llJt Axpd^N X p 



= ReTr H XP exp (-i 3) S% P + S^ 2 P exp [ -i 
+N XP (j: xp exp (-i 3) S$ P - (S 



qA 



<XP ex P 



— 1 



o 



A 

XP 



(41) 



The real part constitutes the desired transport equation, whereas the imaginary part may 
be used to check the consistency of e.g. the approximation for N xp . 

The latter is necessary because, in addition to the approximations made for the retarded 
self energy components, we explicitly need S 12 for the transport equation. Using a similar 
ansatz as for the retarded and advanced self-energy, it is treated as an external field which 
exhibits a jump at t — 0. The values before and after the jump are given by the boundary 
conditions at t — ±oo: 



il2 
J XP 



27rig l0 (9(t)N f T f + 9(-t)N i T i ] 



(42) 



with Nf/i being the final resp. the initial quark distribution function, i.e., the Fermi-Dirac 
distribution at T = Tfu 

It is a rather straightforward task to insert this ansatz into the above equation, after 
some elementary manipulations one obtains the quantum transport equation for the physical 
system described by our aproximations as 

d t N t = 2ng \o(t)T f (N f - N t ) + e(-t)Ti{Ni - N t ) 



-Tig 



Z(t) 



Tf(N f 



-Nj-T^Ni-Nt 
Au 2 t + 7 2 + 7 2 cos(4^t) 



9{t)Z{t) 

- 2ua t sin(4u;ji) 



(43) 



with g as in eq. (pS|). Note, that the real part of the self-energy does not enter the transport 
equation explicitly. The spectral width of the quarks is responsible for the reoccupation in 
the non-equilibrium situation after the jump in temperature. 
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For the interpretation of this result it is useful to study the classical limit of the above 
transport equation given by 



d t N? = 2ng [9(t)T f (N f - Nf) + 9(-t)Ti(Ni - JV t d )J . (44) 

The classical and the quantum transport equation have two things in common. First, if 
the system starts from an equilibrium state, i.e., for t < with N t = iVj (Nf = N), the 
occupation number parameter remains equal to iVj until t = 0. Furthermore, in the limit 
t — > +oo, N t (Nf) approaches Nf. 

Otherwise however, the solutions of the two quations are vastly different. This difference 
between the quantum and the classical transport problem may be characterized by the slope 
of the solution at t — + . The slope of the classical solution at t = 0+ is 2ivgTf(Nf — N), 
while the slope of the full soulution is zero. 

This implies, that the reaction of the quantum system to the jump in temperature does 
not start immediately as in the classical case, but with a certain delay time. In Fig. 3 we 
plotted the classical and the quantum solution for comparison. As may be inferred from 
Fig. 3, the quantum delay time decreases with increasing temeprature. However, the highly 
nonlinear nature of our coupled equations leads to a non-analytical dependence of this delay 
time on the temperature. 

For our model calculation, this behaviour is examined in more detail in Fig. 4, where we 
plotted values of the time r needed to reach 90 % of the final occupation number for quarks 
at rest. The curve parameter is the critical temperature in the range from 150 to 250 MeV, 
marking the boundaries of the range of currently accepted phase transition temperature 
from nuclear to quark matter. Values for r are given only for T > T c . 



VI. CONCLUSION 

To draw conclusions from the present paper, we work backwards starting from Fig. 4 and 
5. The most prominent fact inferred from these figures is the lenghtening of the relaxation 
time for a QGP by a factor due to quantum effects. This result is qualitatively as 

well as quantitatively similar to the results obtained for different systems, cf. [|J. Also, 
the absolute value of the relaxation time is comparable to thermalization times generally 
accepted for quark matter but systematically higher than 1 fm/c. 

In particular: using a phase transition temperature of 160 MeV as obtained in various 
realistic calculations, together with a final temperature only moderately higher than this T c , 
the time r to reach a 90 % thermalized QGP is ~ 3.6 fm/c. In our view, this time scale is 
much too long to sustain the constant "hot glue" scenario. The natural conclusion therefore 
would be, that with a moderately ultrarelativistic heavy-ion collison one might probably see 
a phase transition to quark matter - but most ceratainly not a thermalized plasma phase. 

The quantum effects we find have a straightforward physical interpretation: "collisions" 
need a certain time to build up. However, in our view this straightforward interpretation is 
very misleading. We have shown, that the "quantum delay" is due to the spectral width of 
the quarks in the hot medium. This spectral width is a manifestly non-classical effect, to 
attribute it to "collisions" is sneaking the mechanistic quasi-particle picture back into the 
physical description. 
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That this picture is in profound contradiction to our results is obvious from the quantum 
oscillations in Fig. 2 and 3. Such quantum oscillations must be taken rather serious, because 
in a completely different area of physics they are even established experimentally: They seem 
to play an important role in four wave mixing of ultrashort laser pulses |[23|| . Consequently, 
neither the effective mass, nor the spectral width parameter are monotonous functions. 

To our knowledge there is no principle which would forbid such oscillations also in the 
occupation number parameter N t . Indeed, a close look at the results reveals such oscillations 
in the derivative of our solution, albeit not strong enough in the present model to lead to a 
piecewise diminishing N t . 

With the present paper, we have taken the non-equilibrium character of our quantum 
system more serious than a realistic description of the interaction. In particular, we made 
a physically motivated ansatz for all three independent elements of the 2x2 matrix valued 
self energy function. One might argue, that in a "real" plasma of gluons and, eventually, 
quarks this interaction is much more complicated. This is true in principle, buth it is not 
clear to what consequence: The highly nonlinear nature of the quantum transport problem 
might well lead to an enhancement of the effects we observe. A calculation within another 
scheme, where the real part of the self energy function was treated differently flQ |, indeed 
indicates more quantum oscillations than obtained here. 

A final remark of the present paper concerns the use of old-style classical transport 
equations for strongly interacting systems like colliding nuclei: As we have shown, it is not 
justified a priori, i.e., without checking against a solution of the matrix valued Schwinger- 
Dyson equation. 
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APPENDIX A: SUMMARY OF THE EQUATIONS FOR THE SPECTRAL 

FUNCTION 



In this appendix we present some derivations used above to obtain u t , jt and N t . 
The real and imaginary part of the self-energy can be splitted into a constant part and 
a part being only proportional to a ^-distribution in time. 



ReS 



R/A 
XP 



ReEi + 6(t)(ReE f - ReE;) = ReE; + #(t)ARe£ , 



ilmSj 7 / = T^Txp = T^91o(Ti + (7) - Ti)0(t)) = + AT0(t)) . 

The contributions of the constant parts of the self-energy can be evaluated trivially, when 
inserting the ansatz for the spectral function, eq. ([31]): 

ReSi exp (-i o) S R/A 



ReSj 
ReS, 



E - E'ie 

utlo + P"i + m t 



+ 



^t7o ~P1 -mt 



2co t (E - u t ± v) t ) 2uj t (E + to t ± vyt) 



(Al) 



More care must be applied in the case of the parts of the self-energy proportional to the 
^-distribution. At the same time we evaluate the trace over the Dirac matrices. 



(AReE#(t))exp 



s r/A 



-AReS Tr 

4 



AReS 



9{t) 



E-E'±ie 
1 1 



dE 



6{t) 

m t 

2uj t \E — uJt ± ie E + uo t ± ie 



,e-^ E - E 'H{^t)A{E' ,p,t) 



E-E'±ie 



™ . -2i{E-ui t )t p -2i(E+L0 t )t 

±6( T t)^e^(^ — - ! 

ZUJt \E—UJt±l6 E + UO t ±16 



(A2) 



-Tr 



(Ti^To AT9(t)) exp ( -i $) S R/A 



4 

ingAT 



mg AT Tr 

T0(D 



Tl O 0(t) I dE 



, A(E>,p,t) 
E-E'±ie 
1 1 



- / dE 



f e- 2i ( E - E '»6(Tt)A(E',p,t)- 



E-E'±ie 



2u)t V E — uj t ± ie + o; t ± ie 



e -2i(J5-wt)t 



-2i(£+w t )i 



2uj t \ E — tu t ± ie E + uj t ± ie / . 



+ 



(A3) 
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TABLES 



a S {T c ,T) 


150 


175 


200 


225 


250 


275 


300 T in MeV 


150 


0.313 


0.291 


0.276 


0.262 


0.251 


0.242 


0.234 


175 




0.313 


0.294 


0.279 


0.267 


0.256 


0.248 


200 






0.313 


0.296 


0.282 


0.271 


0.262 


225 








0.313 


0.298 


0.285 


0.275 


250 
T c in MeV 










0.313 


0.299 


0.287 



TABLE I. The coupling g depending on the critical temperature T c and the temperature of the 



quark-gluon plasma T. 



m(T c ,T) 
in GeV 


150 


175 


200 


225 


250 


275 


300 T in MeV 


150 


0.121 


0.137 


0.152 


0.167 


0.181 


0.196 


0.21 


175 




0.142 


0.157 


0.172 


0.187 


0.202 


0.216 


200 






0.162 


0.177 


0.192 


0.207 


0.222 


225 








0.182 


0.197 


0.212 


0.228 


250 
T c in MeV 










0.202 


0.218 


0.233 



TABLE II. The thermal mass m depending on the critical temperature T c and the temperature 



of the quark-gluon plasma T. 
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FIG. 1. tot as a function of time, curve parameter Tj = T c with values 250 MeV (short dashed), 
200 MeV (long dashed) and 150 MeV (continuous). 
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FIG. 3. Classical solution Nf (thin lines) and quantum solution Nt (thick lines) of the transport 
equation, curve parameter Tf = T c with values 250 MeV (short dashed), 200 MeV (long dashed) 
and 150 MeV (continuous). 




0.125 0.15 0.175 0.2 0.225 0.25 0.275 
FIG. 4. Relaxation time r to reach 90 % of the final occupation number for the classical 
(thin lines) and quantum solution (thick lines). Curve parameter is the critical temperature 
T c = 150, 175, . . . 250 MeV, the dots indicate the transition point. 
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